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Abstract 

A variable c(t) has been previously proposed as an alternative to dark 
energy as the cause for the apparent acceleration of the universe found 
with supernovae la on high red-shift galaxies, but without an analytic 
expression for c(t). We have herein derived a variable photon velocity 
from the assumption of anisotropic and homogeneous universe which 
is consistent with being the physical photon velocity c(t) in that it uses 
physical clocks and rulers. The derivation makes c(t) — aE, the cos- 
mic scale factor times the Hubble ratio in normalized variables. But 
a variable c(t) no longer keeps the Lorentz transform invariant. In- 
stead, locally we have a generalized Lorentz transform and generalized 
Minkowski metric. By introducing a generalized time for which the 
differential is c(t)dt, we create a covariant derivative and a Minkowski 
metric in order to use Riemann algebra applied to generalized field 
equations. With this we find that the Hubble ratio has the same de- 
pendence on the cosmic scale factor as it does for a constant photon 
velocity. 

Measurements of the fine structure in clouds of cosmic matter ex- 
perimentally demonstrate that some physical "constants" do change 
with cosmic time. We use dependences on c(t) of various physical 
"constants" which keep a normalized quantum Lagrangian invariant. 
These explain the apparent acceleration of the universe found with 
supernovae la on high red-shift galaxies by assuming a flat universe 
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with no dark (vacuum) energy, and by experimentally determining the 
supernovae la luminosity to be proportional to l/c(t). 

keywords: cosmology theory — distances and redshifts — elementary particles — 
ISM lines and bands — relativity — supernovae la 

1 Introduction 

Measurements of supernovae la at high z [1][2][3] have been interpreted as 
evidence for the existence of dark energy, which has led to its wide spread 
acceptance by the cosmology community, even though there is not a real 
physical understanding of it. There have been suggestions in the literature 
[4] [5] [6] that a variable photon velocity might also explain the experimental 
observations. However, there has not previously been a good analytic expres- 
sion for such a variable photon velocity. This paper reports on a derivation 
of such an expression from the usual assumptions of a homogeneous and 
isotropic universe and the definition and invariance of the space-time line 
element. 

The assumption that the universe is homogeneous and isotropic does not 
automatically lead to a constant photon velocity. In fact, the Robertson- 
Walker (RW) space-time line element [8] [9], which can be derived from the 
assumption of this symmetry, allows for a variable photon velocity c(t) which 
depends on the cosmic time t. Any coordinate system that results from a 
diffeomorphic covariant transform from the RW coordinates will describe 
the same physical system as does the RW coordinates. These constitute a 
geometric manifold with a variable photon velocity. In particular if the co- 
ordinates from one such transform can be read on physical rulers and clocks 
and is variable, then the photon velocity described by its coordinates will 
be physical, and all other diffeomorphic covariant transforms, including the 
original RW coordinates, will also describe physical systems with the same 
physical variable photon velocity. We have discovered a group of coordinate 
systems which have the same variable, observer-independent, photon veloc- 
ity vo(t) at their origin, which, if vo(t) = c(t), we will show have physical 
clocks and rulers. These "physical origin" coordinate systems are a subset 
of all covariant transforms from an RW system (t, x, 0, <f>) constrained by the 
requirement that their spatial points move at a negative velocity when mea- 
sured on the RW radial coordinate and their time T(t, x, 6, <p) be measured 
on clocks attached to the spatial coordinates at a distance R(t, x, 0, <j>) from 
the origin. (In this paper "origin" will always mean spatial origin, where, 
generally, the observer is assumed to be located). 
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But a variable photon velocity has far reaching implications. A large 
number of other physical "constants" depend on it, so they presumably also 
change with cosmic time. The measured change in the fine structure spectra 
of heavy metals at high z demonstrates that such a change does take place [6] . 
In addition a variable c(t) no longer keeps the Lorentz transform invariant, 
so that the equivalency principle and Einstein's field equations must be 
modified. Instead, locally we must have a generalized Lorentz transform 
and generalized Minkowski metric. Following the suggestion of Magueijo 
[7], we introduce a generalized time t for which dt = c(t)dt in order to define 
a covariant derivative and create a Minkowski metric for use in a generalized 
field equation. Of course, the generalized Lorentz transform must reduce to 
the normal Lorentz transform for events whose time separation is small 
compared to cosmic time in order that all the laws of atomic, nuclear, and 
particle physics remain valid. Also, the generalized field equation must 
reduce to Einstein's field equation for the measurements in the Milky Way or 
other close-in galaxies which have been used to confirm Einstein's equation. 

Using Magueijo [7] as a guide, we determine the dependence of other 
physical "constants" on c(t) by making the quantum Lagrangian invariant 
to c(t). This combination explains an apparent accelerated universe which 
we will show agrees with observation of Supernovae la at high z. 

We proceed by using the RW coordinates with c = c(t) applied to a 
radial world line to find the radial 4-velocity and 4-acceleration of a point 
on the frame of the covariantly transformed coordinate system in terms of 
both the RW coordinates and the transformed coordinates themselves (Sect 
2). This will yield partial differential equations (PDEs) for the allowable co- 
variant coordinate transformations from the RW to the transformed system 
(Sect 3). We integrate these PDEs for those transformed coordinates which 
have a diagonal metric (Sect 4). These all have the same variable photon 
velocity vo(t) at the origin, which, if vo(t) = c(t), makes c(t) = aE (Sect 
5). Isotropy requires all of the transforms to be a newly defined generalized 
Lorentz transforms close to their origin (Sect 6). We develop criteria for 
physicality and show that there exists a group of transforms, both diago- 
nal and non-diagonal, which meet these criteria at the origin with the same 
vo(t) as the diagonal coordinates (Sect 7). By introducing the generalized 
time t, a generalized field equation is defined, which locally reduces to the 
generalized Minkowski metric and Lorentz transform (Sect 8). This allows 
us to calculate a(t) and thus c(t) (Sect 9). The dependence of other phys- 
ical "constants" on c(t) is derived from an invariant normalized quantum 
Lagrangian (Sect 10), which then explains the observed apparent universe 
acceleration (Sect 11). We discuss feasibility of earth experiments to mea- 
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sure c(t) in Sect 12. In Sect 13 we point out that c(t) = vo(t) does not 
remove the horizon problem even though c(t) goes to infinity as t goes to 
zero. 

2 The transformation from RW coordinates 

For a homogeneous and isotropic universe, one can derive without the use 
of relativity the RW moving coordinate system whose metric in polar coor- 
dinates (r,0,(f)) [10, page 412] (adapted to c(t)) is 

ds 2 = c(t) 2 dt 2 - a(t) 2 [d X 2 + r 2 duj\ (1) 

where t is the universal cosmic time, x is the co-moving curvilinear radial 
spatial variable (% > 0), and a(t) is the scale factor of the RW metric, 
did 2 = d6 2 + sin 2 6d(j) 2 , and 

{sinx, k = l, 
X, k = 0, (2) 
sinhx, k = -1, 

where k is the RW spatial curvature variable to indicate a closed, flat, or 
open universe, resp.. In order to define a covariant derivative [7], it will be 
convenient to introduce the time related quantity t defined by 

i= [ c(t)dt, (3) 
Jo 

so that the line element becomes 

ds 2 = di 2 -a 2 (d X 2 + r 2 duj 2 ). (4) 

Since this is the exact form of the original derivation of the RW metric, it 
is clear that introducing c(t) does not inviolate the derivation. We will take 
the origin of the transformed radius to be the same as the origin of the RW 
radius: R = at x = 0. We will require that T = t at x = 0, since the time 
on clocks attached to every galactic point is t, including the origin. We will 
also make the magnitude of the photon velocity in the (T, R) system to be 
unity at the origin at the present time To = to- This makes the unit of c(t) 
be the value of what is measured by an observer at the origin at the present 
time (e.g., on earth now). 

We will assume that the angular coordinates 9 S and 4> s of the new trans- 
formed coordinates are given by the angular coordinates of the RW system. 
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We will be initially confining our attention to only radial trajectories, so 
that the angular coordinates are fixed, and T and R will be functions only 
of t and x- Thus 

T = T(t, X ), R = R(t, X ) 

e s = 9, <p s = <p. [ } 

The transform of the differentials is given by 

dT = T t dt + Tydx-, dR = Rtdt + R y dx, 

(6) 

d6 s = d8 = 0, dcj> s = d(j) = 0, 

where the subscripts (other than s) indicate partial derivatives with respect 
to the subscript variable. We will henceforth drop the subscript s. 

We will assume a general form for the metric expressed in transformed 
coordinates with spherical symmetry. [10, page 335]: 

ds 2 = A 2 c(t) 2 dT 2 - B 2 dR 2 - 2Cc(t)dRdT - a 2 r 2 duj 2 . (7) 

where B and C are implicit functions of (T, R) and explicit functions of 
(t, %). We assume that A = 1, which is equivalent to an assumption that 
the clocks on the transformed system are attached to the spatial coordinates. 
Also, since both the RW metric and our transformed metric are spherically 
symmetric, the coefficients of the angular variables are the same. 

We will use as the RW coordinates (t, %, 6, (f>) and as the transformed 
coordinates (T, R,6,(p), where dT = cdT, so that the transformed metric 
becomes 

ds 2 = df 2 - B 2 dR 2 - 2CdRdf - a 2 r 2 du 2 , (8) 
and the transform from (i, x) to (T, R) becomes 

dT = T t dt + cT-ydx, 

X (9) 
dR= ±R t dt + R x dx, 

Since the contravariant velocity vector of a point on the transformed 
system in the transformed coordinates is (1, 0, 0, 0), this will transform from 

the vector in the RW coordinates (7, ^,0,0) the same as dT,dR in eq 
9: 

1 = T tl + lcT xl V, 
= \R tl + \R xl V. ' 
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where V and 7 are defined as 



V = a(t) , 7 = . (11) 

V is related to the ordinary RW radial velocity V of a point on the trans- 
formed frame by 

V = a(t) = c(t)V, (12) 

since V is the RW differential distance a(t)dx divided by the RW differential 
time dt, taken for the co-located point R on the transformed frame. 
Manipulating the second line of eq 10 gives 

V = -^, (13) 

where 

»W = f. (H) 

Thus a(i) has the dimensions of time, and a(t) has the dimensions of length. 
(We use various related non-physical quantities with hats to reduce the 
proliferation of ds). 

If we invert eq 9, put this into the RW metric (eq 4 ), and use eq 10, we 
verify that A is 1, and that the expressions for B 2 and C in the metric (eq 
7) become: 

b 2 = m^ni - (ft) 2 ], 



(15) 

; R x L ' ' aTt 1 ' 

We will consider the contravariant acceleration vector of a test particle 
moving with the same velocity as a point on the transformed frame. This 
covariant derivative of the velocity vector can be expressed in both the 
transformed coordinates 

= = ? + (16) 

and the RW coordinates 

* = € + r?!^ 1 = 7 + = (f ) r + 

& = € + ThAuW + UW) = 7 (^) + 2§ ^ = 7*1 (f ) + > 2 ^, 

(17) 
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where Cf = (^p) > the Ts are the affine connections for their respective 

coordinate systems, U and A are the velocity vector and acceleration vector 
of the test particle, the dot represents differentiation by dt, g = d 2 R/dT 2 = 
d 2 R/dT 2 is the ordinary acceleration for the test particle, and g* = dU T /dT 

It will transform like dT,dR (eq 9): 

A f = i*Tt + cA*T y , 

X (18) 
A R = lAtRt + A*R X . 

We can obtain a relation between g, g*, and C by setting to zero the scalar 
product of the covariant velocity vector and the contravariant acceleration 
vector in transformed coordinates: 

CC ~ C ~ 

g^A" =g* + - C{g + ^f^) = (19) 

so that 

g* = Cg. (20) 

If we insert the value of g* from eq 20 and C from eq 15 into eqs 16 and 
17 and rewrite the two equations of 18 in terms of their factors, using eq 13, 
we obtain 

C- R ^ ' 

9+ = -f F ^ 



where the factors are defined as 

(22) 



3 The partial differential equations for a diagonal 
transform 

We are interested in seeing whether there is a diagonal solution (C = 0) of eq 
21, and if there is, what are the PDEs that determine its dependent variables. 
(In this paper we will use the adjective "diagonal" to describe any quantity 
associated with a metric which is diagonal when expressed in transformed 
coordinates, such as diagonal solution, diagonal metric, diagonal transform, 
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diagonal coordinates, etc.). We can see from eq 15 that when F\ = 0, then 
(7 = 0, so there is diagonal solution of eq 21 for iq = 0, g = -^Fi- 

We will now find a differential equation with only V as the dependent 
variable. Transposing F\ = yields another expression for V: 

v = ~W t - (23) 

Putting eq 23 into eq 10 produces a simple relation for T t : 

T t = 7 . (24) 
Now we can write a formula for T, using eq 23 and eq 24: 

T = t+ T x d X t = t+ (-ajV)dxt, (25) 
Jo Jo 

where the symbol d\t signifies that the integral is to be carried out with 
respect to \ with t held constant, and we have used the boundary condition 
that at % = 0, T = t. Eq 25 can be partially differentiated with respect to t 
(giving 7) and then with respect to x an d use eq 11 to obtain a PDE for V: 



R 



R 



We will find these diagonal transforms very useful since they represent 
the simplest covariant transform. We will be able to find a general solution of 
the PDE's for them, including a formula for calculating the variable photon 
velocity at the origin vo(t). This is very convenient, since we will find that 
there is a physical origin group of coordinates, diagonal and non-diagonal, 
all of which have this same vo(t) (see Sect 7). 



4 General solution of the partial differential equa- 
tions for a diagonal metric 

Thus, we shall now try to find expressions for the transformed coordinates 
with a diagonal metric. The PDEs we need to solve for T(t,x), R{t,x), 
and V(t, x) are the three equations 13, 23, and 26 subject to the boundary 
conditions at x = of R = 0, T = t, V = 0, and the photon velocity equal 
to one at t = to, x = 0. Eq 26 can be rewritten as 

9Vr = _9or 
V(l-V 2 ) a ' 
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where the subscript on the partial differential indicates the variable to be 
held constant. This can be integrated with an integration constant ln/c 
(k > 0). Since the integration is done at constant R, then k = k(R), and 
inversely, R = R(k). Integrating eq 27, we get 

* = (28) 
va + k 

where we have selected the minus sign for the square root in accordance 
with our assumption that the transformed frame have a negative velocity 
when measured in RW coordinates. 

At this point, R is an unknown function of k. The various possible 
coordinate systems which solve our PDEs are characterized, in large part, 
by the function R(k). But for all, in order for V to vanish when R = 0, k 
must also; so always 

k(0) = 0. (29) 

One observation we can make without any more work is that, as long 
as k(R) remains finite, V goes to —1, and V goes to — c(i) (see eq 12), for 
a(t) = 0, i.e. for t = 0, the horizon. This is different from any of the three 
definitions of distance assuming a constant photon velocity at the origin 
treated by Fletcher (including the proper distance) [11], for all of which V 
became —1 at a finite t instead of at t = 0. 

Let us now look at lines of constant k(R), i.e. constant R, in t,x space. 
Eq 11 can be integrated for \ with use of eq 28 at constant k to give the 
following general solution, for an arbitrary function f(n) and an arbitrary 
upper limit t*: 

x(t,n) = f( K )+ T — ffi = . (30) 

f(n) is the value of x(t, n) for t = t*. For an open universe, we expect that 
if R is kept constant the galactic point \ that will be passing any given R 
will eventually approach zero as RW time t approaches infinity. Thus for 
this case we can set /(ft) = at t* = oo. 

For a closed universe, a will be periodic, so that V goes to zero at some 
point, which we will take to be t* , after which we may choose a positive sign 
for V to represent a collapsing universe. (Actually, for this derived c(t), if 
we solve the field equation in Sect 8, we will find t* = oo, so the universe 
never collapses). In order for x to remain finite at that point and still give 
physical coordinates at the origin (see Sect 7.1), we will make 

f{K) = M f ^/tu r (31) 

Jo a{s)^/a z (s) + k z 
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Here, the constant M will be a function of the universe energy density, 
so that it will become zero at the transition between a closed and a flat 
universe. 

At this point, we have obtained V = V*(t, re) from eq 28 and have also 
obtained the function %(i, re). We can in principle invert eq 30 to obtain 
re in terms of t and x : K = K(t,x)- This gives us the velocity function 
V(t,x) = V*(t,K(t,x)- If the function R(n) were known, we would then 
also have R(t,x) = R(K(t,x))- 

T(t, x) can be found by noting from eqs 23, 24, and 28 that 



T x = -VaT t = -Va/\] 1 - V 2 



K. 



(32) 



By substituting eq 32 into eq 25, and integrating over re instead of x by 
dividing the integrand by the partial of eq 30 with respect to re, we find an 
expression for T(t, x)- For an open universe it can be written as 



T(t, X )=t + 



s: 



i - 



a(s) 



A/a 2 (s) + k 2 



ds K , 



(33) 



where re is put equal to K(t, %) after integration at constant re in order to 
get T(t,x)- For a closed universe we can use eq 31 to get /(re) and thus 
T(t, X ): 



T(t, X )=t + M 



a(s) 



1/ a 2 (s) + re 2 



ds K + 



a(s) 



^/a 2 (s) + K 2 



ds K . 



(34) 



This completes the solution. These are all the possible solutions for the 
diagonal case. We note that these are in reality all the same system except 
for different calibrations of the R axis with R = i?(re). 



5 Photon velocity close to the origin for a diagonal 
metric 

We would like to find an expression for photon velocity in the transformed 
diagonal coordinate system at the origin, that is as re approaches zero by eq 
29. Using eqs 13, 23, and 24 inserted into 15, the metric (eq 7) becomes 

ds 2 = c 2 dT 2 _f^L\ 2 dR 2_ a \1 dbJ 2 . (35) 
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Although B(t, x) has been determined from radial trajectories, when 7 and 
R x are expressed in terms of t and %, the metric of eq 35 is valid for all 
trajectories, since for spherical symmetry, B cannot depend on the angular 
coordinates. The radial photon velocity is given by dR/dT for ds = 0: 



R x 
7a 



(36) 



(For C = we write the positive value of the photon velocity, but understand 
that in all cases for a diagonal metric there are both positive and negative 
values). 

From eqs 36 , we get 



v P (t,x) 



1 dR{n) 



"fa 



dn 



dK 



(37) 



K x can be written in an inverted form by taking the derivative of eq 30 with 
respect to k: 



( dK 



a(s)ds K 



t (a 2 (s) + K 2 ) 3 /2' 



(38) 



We will let vo(t) be the photon velocity for a diagonal metric at the 
origin. To find it, we set k = in eq 38, 7 = 1, and use eq 37 to give 



v (t) 



ok'(0) 



ds 



(39) 



k'(0) is a constant to be adjusted to make vo(to) = 1. It should be noted 
that the value of vo(t) is independent of the functional form of k(R), and is 
therefor the same for all diagonal transforms with the same a(t). 

We will be showing (sect 7) that it is consistent to have the physical 
photon velocity c(t) equal vo(t). Under this assumption, eq 39 with a change 
in the integration variable from t to a(t) eq 39 becomes 



- = *'(0) 
a 



f'(o) + f 

J a 



c(a)da 
a 2 a 



(40) 



remembering that the dot indicates differentiation by t. By differentiation 
of both sides of eq 40 by a, we can obtain 



c(t) 



1 



a = aE, 



(41) 
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where a is the normalized scale factor 



a = a/ao, (42) 
and E is the normalized Hubble ratio H(t) 

H 1 a 

v^jr^jr-- (43) 

(The subscript denotes the value at t = to, the present time; k'(0) = clqHo). 
From this relation, the generalized field equation (sect 8) will enable us to 
evaluate a and E and thus c(t). 

6 Transformed metric near the origin 

Since V goes to zero as R goes to zero, and we are assuming isotropy of the 
universe, the photon velocities in all directions become the same, with C 
going to zero at the origin for all the transformed coordinates. 

The metric for all the transformed coordinates for small R can thus be 
written 

ds 2 = c 2 dT 2 ^— [dR 2 + R 2 dui 2 ] , (44) 

Vp(t, 0) 

where v p (t, 0) is the variable photon velocity at the origin, which is given by 
vo(t) in eq 39 for diagonal metrics, but may possibly be different for some 
non-diagonal metrics. 

In this limit, since V(t, 0) = and C(t, 0) = = j£(t, 0) and T t (t, 0) = 1 
(eq 10), then from, eq 15 

1 -^fl) = AA, (45) 



v p (t,0) c y ' R x (t,0) R 

where the last equality has used the first terms of the series expansion of R 
around x = [Rt(t, 0) = 0]. Therefor, all transformed coordinates near the 
origin become 

T = t 

R = v p (t,0)a(t) X . K } 

Note that the value of R at T = to is given near the origin by a(to)x 
(since v p (to, 0) = 1, a(to) = a(to)), the expression for which all measurements 
of distance reduce for small x [Weinberg 1972 page 424]. 
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If v p (t, 0) = c(t), we can call eq 44 a generalized Minkowski metric 
because for a range of times, ST « t, the metric becomes the Minkowski 
metric. It can be written with the generalized time T: 

ds 2 = df 2 - [dR 2 + R 2 dw 2 }, (47) 

It is the metric for a generalized Lorentz transformation between a frame 
(dT', dR') moving at a radial velocity v relative to a second frame (dT, dR): 

dT' = i(t)(dT - -fadR), 

VA ^ h (48) 
dR' = i{t){-vdT + dR), 



where j'(t) = 1/y 1 — v 2 /c(t) 2 . Of course, because of spherical symmetry, 
this generalized Lorentz transform applies to a frame moving in any direction 
for differential distances measured in the direction of the velocity v. 



7 The Physicality Criteria 

We would like to find those transformed coordinates which are physical since 
this is the type of coordinate system with which our physical laws are for- 
mulated. The cosmological principle gives us a unique opportunity to do 
this, because it says that the clocks and rulers on every galactic point, when 
properly calibrated and synchronized, are the same. So, to be physical we 
need the transform to use the same clocks and rulers and to have the same 
photon velocity as the co-located galactic point. This latter requirement 
arises because at every galactic point we have the generalized Minkowski 
metric for which the photon velocity is independent of the observer, includ- 
ing co-located observers on the transformed frame moving with a velocity 
V. For criteria of sameness, we apply a version of the equivalency principle 
that the laws of physics are the same in every physical coordinate system. 
We assume measurements in every physical coordinate system made on a 
clock stationary in that system are all the same, and measurements on a 
ruler for constant time in that system are all the same, for a universal line 
increment ds. 

We are already assured that the clocks which measure T are the same 
as the clocks at the origin by making the metric coefficient of dT be unity, 
since then c(t)dT(t, R) = ds at constant R, and c(t)dt(t, x) = ds at constant 
X, so that the attached clocks on the transformed coordinates are the same 
as the clocks at every galactic point. From the co-located galactic point 
having earth rulers, we get a physicality requirement on rulers in the radial 
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direction, by making the transformed radial distance 8Rt at a constant time 
T on local transformed clocks be identical to the radial distance a(t)dxt ( e< l 
14) on a co-located galactic point taken at a constant time t on the local 
galactic clocks. This must be true for a given differential of proper distance 
ds, so the radial rulers will be the same if 



dR\ (d X 
= a(t) 



8S ) r 



ds 



(49) 



By putting the left side of eq 49 into eq 7 and the right side into eq 1, we 
obtain the radial ruler requirement B = 1, or if B = B/c, then 



B(t,R) = 



C (ty 



(50) 



This must be independent of R since the clocks, rulers, and photon ve- 
locities on co- located galaxies are all the same at the same cosmic time t, 
independent of R. 

Since B and the photon velocity are generally a function of R, none 
of our transformed coordinates meet the physicality requirement for all R 
except for the zero density universe (see Sect 9). Another exception is at 
the origin, where we will find a subset of the transformed coordinates that 
can meet the physicality criteria for a small range of R. 



7.1 Physical origin for a diagonal metric 

We can make the transformed coordinates physical for a limited range of R 
near the origin of a diagonal metric by requiring that the partial derivative 
of B(t, R) re R vanish at the origin. This satisfies both the ruler requirement 
and the constant photon velocity requirement, since v p = 1/B for a diagonal 
metric. Eqs 35 and 38 show that 



B(t,R) 



07 
R v 



07 



dt' 



I V R (t',R) 
Jt a{t 



(51) 



By partial differentiation, we obtain 



Br = ^— + a7 

-fl v 



fRR- V RR (t',R)—Z. 
Jt a(t'. 



(52) 



Since 7j?(t, 0) = 0, V RR (t,0) = will make /r_r(0) (eq 31) and the square 
bracket of eq 52 and thus B R (t,0) vanish at the origin. If we partially 
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differentiate eq 28 twice, we see that Vrr(£, 0) = requires that krr(0) = 0. 
Thus, those diagonal transforms which have R linear in n for small R have 
physical clocks and rulers which measure the physical photon velocity near 
the origin if 

c(t) = l/B(t,0) = v (t) (53) 

given by eq 39. This makes c(t) = aE (eq 41). 

These criteria are consistent with the more general criteria developed by 
Bernal et al[15]. In their notation for x\ and x 2 translated to the notation 
in this paper, if we assume the RW metric has the length measure dx\ = 
a(t)dx and the transformed length measure is dx2 = dR, then at the origin 
d 2 x 1 = a(t)xR = a(t)/R x = c(i)/v (t), and dix 2 = R x /a(t) = v (t)/c(i) . 
Since these two partial derivatives are equal at the origin when vo(t) = c(t), 
the two systems meet the Bernal criterion for using the same rulers there. 
A similar relation applies to clocks in the two systems, which are the same 
in both systems for all R. 

If vo(t) is to be considered physical, every other transform that meets 
these same physicality requirements at its origin must have the same vo(t) = 
c(t). There are, of course, many coordinate systems, including transforms 
from the RW system, which do not have physical coordinates anywhere. In 
order to find another transform which is physical at the origin, we will try 
to find a physical origin for a non-diagonal transform from eq 1 in an open 
universe. 



7.2 Physical origin for a non-diagonal metric in an open uni- 
verse 

For a non-diagonal metric, the criterion at the origin for constant rulers 
requires that Br = 0. We will introduce the symbol 

T = - ^r, (54) 
aT t 

so that B in eq 15, usng T t from eq 10, becomes 

= 1TT^?- (55) 

R x 1 — VT 

Eq 11 is still valid for non-diagonal metrics, so for an open universe 

X 



r°° - Fit' 

/ V(t',R)^L (56) 
Jt a{f) 
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By partial differentiation of B (eq 55), noting that V(t, 0) = and since 
C(t, 0) = = T(t, 0) = 0, we get 

B fl (t,0) = a(t)(^^) (t > 0) = -a(t)jf t ^(f.O)^ (58) 

Thus, one general requirement for making an open metric, either diagonal 
or non-diagonal, have radial rulers equal galactic rulers near the origin is to 
set Vnii(t, 0) = 0; i.e., have V oc R to first order in R. 

The constant photon velocity imposes the additional condition Cn(t, 0) = 
0, which is equivalent to Tr(£, 0) = Vjj(t, 0) (eq 15). Even though this re- 
quirement makes the metric diagonal to first order in R, there is no a priori 
reason that it have the same vo(t) as the diagonal metric, although we will 
find that in fact it has. 

We need first to find an expression for the transformed time for a non- 
diagonal metric. From eq 10 we obtain 

dT\ V 1 

m) R = Tt + T *W) = WJi)- (59) 

Since eq 59 is valid for both diagonal and non-diagonal coordinates, the 
integration constant is obtained by requiring T be the diagonal expression 
(eq 34) when 7(5, R) equals \/l + K 2 /a 2 . Thus the integral of eq 59 becomes 

ft o./ ,-oo 

T ee G(t, R) = ^ + ^ [1 - l/ 7 (t', R)] dt' R . (60) 

Note that G(t, 0) = t since j(t, 0) = 1. If we partially differentiate eq 60, we 
get Gji(t, 0) = and partially differentiate a second time we obtain 

Vl(t',0)dt' R . (61) 

To find an expression for T, we note from eq 10 that T t = 1/7(1 — VT) 
and by the chain rule that T x = GrR x . Inserting these into eq 54, after 
some manipulation we get: 

~ : (62) 



-V + a/ 1 R x G R 
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The partial derivative of eq 62 at the origin then becomes 



T "" |01 " iW/«»P) (63) 

We can obtain a solution by making the integrands of G RR (eq 61) and 
l/R x (eq 57) be proportional with a proportionality constant of k'(0). This 
makes 

Vh(*,0) = -^M (64) 
a(t) 

and G RR (t, 0)R x (t, 0) = k'(0). Then from eq 63 we get 

r R (t,0) = -^ = V R (t,0), (65) 

showing that eq 64 is a solution. Thus the physical origin non-diagonal 
photon velocity in an open universe becomes the same vo(t) as the diagonal 
coordinates (use eq 64 in eqs 57 and 47), which it must if it is to be considered 
equal to c(t). It seems likely that this is also true for a closed universe. 

It should be observed that up to this point in the paper, we have not 
used the field equation. The RW coordinates, the transform from them, the 
variable photon velocity at the origin, and the physicality criteria depend 
only on the assumed symmetry of homogeneity and isotropy of the universe 
and the invariance of the space-time line element ds 2 . 



7.3 Non physicality of the RW proper coordinates 

The customary concept of radial cosmic distance is the RW proper distance, 
defined as a(t)x, which for constant x g° es to zero as t and a(t) go to 
zero. If we use this distance as part of the RW proper coordinate system, 
R* = a(t)x, T* = t, the metric of eq 1 using eq 14 becomes 

_ f-R^dJ 1 . (66) 
X 2 

This system is the reverse of our transformed coordinates in that the ruler 
physicality requirement is met for all distances but the clock physicality 
requirement is not met at any distance, except at the origin. The coefficient 
A* 2 of dT* 2 is given by 

A* 2 = c 2 - ^(^) 2 (67) 



ds 2 = c 2 dT 



*2 



dR * _ £L^r* 

a dt 
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At the origin, A* 2 = c 2 to first order in R* , so both the clocks and the rulers 
meet the physicality requirement at the origin. 
The radial proper coordinate photon velocity 

/?* dn 

v ;(t,R*) = ±c(t) + —-. (68) 
' a at 

At the origin v*(t, 0) = c(t) and the metric becomes the generalized Minkowski 
(since near R* = 0,x = r )- However, v*(t,0) is not constant to first order 
in R. Thus v*(t,0) does not meet our physicality requirement for constant 
photon velocity. The coordinate photon velocity becomes the physical pho- 
ton velocity at a point on the R* axis, the origin, because of its intrinsic 
composition, not because it is derived by its equivalence to a constant galac- 
tic photon velocity. c(t) is not determined by v*(t,0), and may assume any 
function of t. 



8 Revised field equation for a variable photon ve- 
locity. 

If a variable photon velocity c(t) is a physical reality, we need to revise 
the Equivalency Principle and the field equation to reflect the expecta- 
tion that all non gravitational physical measurements, anywhere and any 
time, obey the generalized Lorentz transforms (eq 48) with the generalized 
Minkowski metric (eq 47) rather than the usual Lorentz transforms and the 
usual Minkowski metric. The use of t or T in the field equation accomplishes 
this: 

87T"(jr 

G^ + Ag^ = -^T^, (69) 

where G M „ is the Einstein contracted curvature tensor determined from the 
coefficients g^ v multiplying the differentials in the metric of eq 4, or of eq 8, 
A is the cosmological "constant" possibly representing some kind of vacuum 
energy density, G is the gravitational "constant" and T^ u is the stress-energy 
tensor. Locally, this will give the generalized Minkowski metric (eq 47) and 
generalized Lorentz transform (eq 48) for dT and dR close to the origin, as 
well as for dt and dd\- 

It will be convenient to use the variables t, r, 8, <f> in the generalized RW 
metric (eq 4): 

ds 2 = d p _ a 2 (^-^ + r 2 d6 2 + r 2 cosW). (70) 

1 KT 
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so the g^ v of eq 69 are the coefficients of eq 70, exactly the same as the usual 
RW metric with constant c = 1 in the variables t, r, 6, <fi. We will make the 
usual assumption that the universe is an ideal fluid with an energy density 
of pc 2 and pressure p, so that we can write the two significant field equations 
[12, page 729] for a(t) as 

3d 2 3k A 8vrG 

— + — -A=— r p, 71 

and 

a a? k SttG , . 

+2- + — + — - A = -p, 72 

a a A a z c 4 

where the dots represent derivatives with respect to t. 

Now, the first equation can be multiplied by a 3 /3, differentiated, and 
subtracted from da 2 times the second to give 

d Gp 3 3G . 2 a 3 - 

Jt^ a) = -^ aap -^ A > (73) 

where we have allowed for the possibility that G and A may be functions of 
c(t). For small p and A, 

Gpa 3 /c 2 = GWH/c 2 ,, (74) 

a constant. If the energy density consists of n particles per unit volume 
of mass m, so p = nm, then the conservation of particles requires na 3 be 
constant. This makes 

Gm 

— tt- = constant. (75) 
c z 

Even when p and A are not negligible as in the early universe, it is reason- 
able to assume that matter particles are conserved well into the relativistic 
particle era and that the relation of constants in eq 75 will continue, so 
that for matter density eq 74 remains valid. But for radiation, we assume 
kT oc 1/a as it is for c = 1 so that radiation energy density is proportional 
to 1/a 4 (see sect 11). 

9 Calculation of a(t) and c(t). 

Following Peebles [13, page 312], we define 



n = -TTTT^T ( 76 ) 
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and 



and 



n r = * (77) 

n Q a Q 

n* - 3^. (78) 
If fcT oc 1/a as it is for c = 1, then we can introduce a radiation term: 

n * = |^^o(T ao) 4 (79) 

where So is the Stefan-Boltzmann's constant divided by c and multiplied by 
1.68 to account for neutrino energy density [13, page 164]. Eq 97 will show 
that B is independent of c(i) so that £Ir 1.7x10 , a constant. Radiation 
makes an important contribution to a only when a < f2#/f2. The £ls are 
defined so that 

n R + n + n r + n A = i. (so) 

Then the normalized Hubble ratio E in eq 43 is determined by eq 71: 



V a 4 cr a z 

which allows us to evaluate c(t) = aE in eq 41. At t = to'- ao = 1, 
-E(ao) = 1) an d Co = 1. 

The physical time t becomes 

where cq = 1 is written explicitly for clarity. For = 1: c = a -1 / 2 = 
(to A) 1/4 , and c #o* = a 2 / 2 , whence c H t = 1/2, a/a = (i/i ) 3/4 - 

It is interesting to consider the limiting case of a zero density universe 
(Q = 0, Q r = Eq 41 makes c = l,i?o«o = 1 and eq 82 makes 

a = i?o^) a = a = t. Eq 30 gives x = csch" 1 (t / k) , or if(f, %) = tsinh\- We 
can then find from eqs 28 and 11 that V(t, x) = ~ tanh \ and 7 = cosh \ so 
that the photon velocity for all distances by using eq 37: 

dR K y dR , . 

v p = * = — 83 

dn jt dn 

Thus the physicality condition is met for all R with R = K and v p = v = 
c(t) = 1, so that the complete transform with eq 34 becomes 

R = tsinhx, 

T = tcosh X - { ' 
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This solution has been known ever since Robertson [14] showed that this 
transformation from the RW co-moving coordinates at zero density obeyed 
the Minkowski metric. What is new is that this solution was derived from 
the equations we obtained for our covariant diagonal transforms and our 
physicality requirement applied for all distances. 

For experiments attempting to measure the variation of the photon ve- 
locity at the origin at the present time, the derivative of c(t) (eq 41) will be 
more useful: 



1 



1 dc 

c dt 



= i-^n-n r = ~ + n A . (85) 
t=t 2 2 



Notice that this fraction is negative when matter dominates, and increases in 
magnitude from zero at zero density to —1/2 at the critical universe density. 
A vacuum energy density opposes the gravitational effect of matter, and as 
it increases, it causes the variation to go to zero before reversing it from a 
decreasing function of time to an increasing function. 



10 Other physical constants. 

First, we should recognize that the spectra of atoms and molecules can 
depend on c(t). Thus, the fine structure constant a j in SI units [16] is 

= sss? <86) 

and the Rydberg constant expressed as a frequency is 

2 "T-eC 2 e 4 m e 

Roo = a f^n = 4(^nr (87) 

which both can depend on c(t) since if c(t) depends on time t, we must allow 
for the possibility that eo will also depend on it because c(t) 2 = l/eo(t)fio(t). 
[Note that the 4neo is often omitted in the fine structure constant since it is 
unity in Gaussian coordinates, but it is essential here if we are to consider 
a variable c(i)]. 

Recently, observers have measured a shift in the fine structure in ab- 
sorption spectra of matter between distant quasars and us [6]. For a range 
of cosmic times, they attribute this to a change in aj for which they found 
Aa f/ctf = -0.72 ± O.I82IO" 5 over the red shift range 0.5 < z < 3.5. 
This experimentally demonstrates that some physical "constants," must be 
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changing with cosmic time. However, the magnitude is far too small to be 
caused solely by the c(t) calculated here. This implies that if c(t) = vo(t), 
some other physical "constants" must vary as well, such as those found in 
eqs 69, 86 and 87. 

A variable physical photon velocity has been previously investigated for 
general consistency with all our physical laws [7] with suggestions for how 
this might explain many cosmological puzzles [4] [5] . We are here concerned 
with a derived variable physical photon velocity and not with assumed vari- 
ations to explain these puzzles. We here propose a possible variation of 
other physical "constants." Following the suggestions of Magueijo [7], we 
can obtain two relations: (1) To make the interactive term of a normalized 
quantum Lagrangian independent of c{t) for particles with a charge gi in a 
field (gi = e for electromagnetic field), we make gi/fic constant; (2) To make 
the matter term independent of c(i), we make m^c 2 /he constant, where rrii 
is the mass of a particle of i matter. We can put these relations together 
with eq 75 and the experimental observation that aj is almost constant to 
obtain a parameterized dependence on c q : 

c 4 1 

gi,hc,miC 2 ,—,e , ~ oc c q . (88) 

G hoc 2 

It is possible that Magueijo's minimal coupling requirement q = applies, 
for which the rest energy of individual particles is independent of c(t). This 
will make the quantities in eq 88 all constant. Unless otherwise stated in 
this paper, we will take q = 0. If we apply eq 88 to an emitted spectral line 
Vnm (eq 87), we get 

fnm OC 

i?oo oc c(t), (89) 

independent of the parameter q. The small change in the separation of spec- 
tral frequencies determining an apparent shift in ctf with c(t) in intergalactic 
gas could then be considered a secondary effect as the ratio of aj to 
decreases with c(t). 

If we review the evidence for unchanging constants [12, sect 38.6], the 
arguments for no change with spatial location is not violated by the present 
determination that they change only with cosmic time t. All the arguments 
for no change with time seem to be related only to the fine structure constant 
a/, which we are also assuming doesn't change with time. 

The change in u nm (c) will change the red shift formula [10, page 416, 
using a]. Although the elapsed time At e and the photon period are each 
increased by ao/a e by the expansion of the universe, the frequency starts 
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out at c(t e ) times the normal frequency at the origin so that 



1 + ^ = 7^- = ^ = -, 00) 

c[t e )a e a e a e 

where t e is the cosmic time of emission. For = 1 this will make c(t) = 



11 Distance determination from supernovae la 

The computed distance di [10, page 421] from a standard candle is derived 
from the received power/area being reduced from the emitted power/area. 
Because of the expansion of the universe, the photon freqency is reduced 
by a e / a o an d the time beween the emission of photons is increased by the 
same ratio. Distance, including space curvature, further reduces the power 
per unit area by 47r(ao?" e ) 2 . Planck's constant increases the energy of each 
photon by c e to give a net ratio of received power/area to emitted luminosity 
of (a e /ao) 2 c(t e )/(47raQrg). This results in an apparent distance given by 



d L = 



o' e 



,L 



(^) 



Ce a e 



V2 =ao(1 + z)c l/2 re( ^)l/2 



(91) 



where we have allowed for the possibility that the emitted luminosity L e of 
the standard candle might change compared to its value L s determined in 
local galaxies at the present time. 

If we expand c e and r e in powers of z, we get 



r e = sink 



aE = l + f(ft-2f) A ) + O0z 3 ) 

(i + \n - n A ) + o(z 3 ) 



r*o cdt 


1 


Jt e a 





2 



(92) 



In order to allow for variation in L e , we will parameterize it by making it 
proportional to 1 + z) n : 

L e /L s = (1 + z) n . (93) 
di thus becomes to second order in z 



dL = 7T 



z + y(-n + ^ + ^ r + ^a) + 0(z 3 ) 



(94) 



High z measurements of supernovae Type la [1] [2] [3] indicate an appar- 
ent accelerating universe. They have shown that the best fit to their data 
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using a constant photon velocity is obtained for a fiat universe (fi r = 0) by 
$7 = .3 ± .15, = .7 =f .15(68.3% confidence). For a constant c = 1 the 
power series expansion of d* L (z) becomes 



This gives a coefficient of z 2 of 0.78 ± .11. Using eq 94 with c(t), we get 
the z 2 coefficient to be 0.75 with Q = l,n = —1/2 (Q\ = Q r = 0). That 
is, it appears that we can explain the observed anomalous acceleration by 
the variable photon velocity calculated for a flat universe without requiring 
the existence of vacuum energy. However, it requires the maximum emission 
of the supernovae Type la to be proportional to 1/(1 + z) 1 / 2 = l/c(t). It 
should not be surprising that the peak luminosity of the SNIa varies with 
c(t). Thus, for a hydrogen-like atom, the emitted power also varies with 
c(t), but in this instance being proportional to c(t). For the SNIa it will 
take a theoretical analysis of the physics using c(t) to see whether it agrees 
with the experimentally determined l/c(t). 

Recently, the Supernovae la have been extended to even higher z[17]. For 
high z, the observed spectrum shifts to longer wavelengths by a factor of 1+z 
compared to the spectrum at low z (by definition of z). As experiments reach 
to these higher values of z, comparison of the derived di with experiments 
is a little subtle. Thus, the correction to the rest frame of the SNIa means 
that the measured times should be divided by ao/a e = E instead of 1 + z 
used in the paper. For 0, = 1 this correction is E = (1 + z) 3 / 2 . Also, we 
should use the full expression for di in eq 91 instead of a power series in z 
(eq 94). For f2 = 1 and n = —1/2 this becomes 



In addition, the spectrum in the rest frame of SNIa increases in frequency 
by c(t e ) compared to the low z spectrum, so that the wavelengths in the 
rest frame are the observed wavelengths divided by E/c(t e ) = 1 + z, just as 
assumed in the Tonry paper [17]. Preliminary reworking their measurements 
with these adjustments indicate that they agree with eq 96 within the quoted 
uncertainties. 

The success of this comparison suggests future work with c(t), especially 
where vacuum energy has been previously invoked, e.g., CMB, gravitational 
lensing, and dynamical estimates of galactic cluster masses. 

Another conclusion we can draw from q = is confirmation that kT for 
decoupled radiation varies inversely with a(t) (and not with a(t)). Thus, if 



i- z+ Z -(l-Ul + Sl K )+0{z*) . 



(95) 



H d L = 2(1 + z)(-l + s/l+~z) 



(96) 
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the radiation starts out as a Planck distribution of photons with frequency, 
each frequency will be decreased by 1/d. so that the distribution is main- 
tained as the universe expands (as is observed) with the maximum of the 
distribution occuring at then kT oc hv max oc 1/ca = I /a. There is 

no separate dependence on c(t) of k or T; the determination of the units 
of the temperature scale determines k. It should be observed that when we 
consider radiation in the early universe that the value of Bq in eq 79 is 



'60 {hcf 



B o = 1.68-t-^, (97) 



a constant for q = 0, so Q,r is also constant, as was assumed in the formula 
for E(a)(eq 81). 

It would also be a test of the theory if the secondary effect of the fine 
structure spectrum in the heavy metal atoms at high z could be calculated 
and compared to observation. Another challenge is to calculate the age of 
old globular clusters using c(t) with fi = 1 to see whether it is shortened 
enough to be significantly less than the 16% shorter age of the universe, 
1/(2H ) instead of 2/(3H ). 



12 Feasibility of confirmation by earth based ex- 
periments. 

A natural question is whether the predicted variability of a physical c(t) is, in 
principle at least, subject to confirmation by earth based experiments. For 
a critical universe density with Hq = 50 km/s/Mpc and no vacuum energy, 
so that c(t) = (to/t) 1 / 4 , we would have a fractional decrease in the measured 
photon velocity of At/4t = .8a;10~ 14 in 10 4 sees. We note, however, that if 
we have vacuum energy density, that this will affect the photon velocity in 
the opposite way, so it could cause an increase in the photon velocity (Eq 
85), rather than a decrease. 

Recent developments project that fractional frequency measurement ac- 
curacies of spectral emissions should achieve 10~ 17 over a period of 10 sees 
in the near future [Berquist et al 2001; Kleppner 2001]], which is a thou- 
sand times more accurate than the predicted cosmic variation of c(t). If our 
proposed variation of physical "constants" described in the previous section 
is valid, then v oc c(t) and A = c{t)/v = constant, so that measuring the 
change in v would be equivalent to measuring the change in c(t). 
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13 Horizon of CMB 



We would like to see if the increasing c(t) when t goes to zero allows com- 
munication between CMB coming to us from different directions. The CMB 
is characterized by the microwave frequency reaching us has a red shift of 
about 1 + z = l/a e = 1000. In the time from the beginning of the universe, 
photons could transverse laterally (constant x) a distance A of 



<e / x , I [ ae da 



Along a transverse angle 6 e in the observer's coordinate system, this distance 
is given by A = a e r(xe)0 e for small angles, where r(x) is given by eq 2, and 
Xe is given by 

[ t0 cdt - 1 f 1 da 
Xe Jt e a a H J ae a 2 E 

For a flat universe, where r(x) = X> t ne angle of the lateral horizon is given 



by 



f«e da 
JO aE 

a f 1 -4?L 



(100) 



1/2 

For Q = 1, a e » Qr and 9 e = a e /3 = .01 radians. Only for angles 
smaller than this are the CMB ripples coordinated (except by inflation). 
Inflation is still needed to establish communications in the early universe. 

A rather surprising result is that this horizon is exactly the same as for 
a constant c = 1. In that case, dt = cdt becomes dt and a(t) becomes a(t) 
so the equation eq 100 is the same. 



14 Conclusions 

Using the assumption of a homogeneous and isotropic universe and the in- 
variance of the space-time line element, we have derived a group of "phys- 
ical origin" coordinates, covariantly transformed from the RW coordinates, 
which at their origin are consistent with being physical with a general- 
ized Minkowski metric and an observer-independent variable photon velocity 
vo(t), the same for the whole group, which includes metrics that are both 
diagonal and non-diagonal. By consistent, we mean that if their photon 
velocity equals the physical photon velocity and is constant to first order in 
R at their origin, then they use the same rulers and clocks that we have at 
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every galactic point. We have here explored the possibility of their being 
equal, which determines c(t) to equal the cosmic scale factor times the Hub- 
ble ratio: c{t) = aE in normalized variables. We have also shown that the 
coordinate photon velocity of the RW proper coordinates is not constant to 
first order in R at the origin and does not determine the physical photon 
velocity c(t). 

Assuming we know this physical photon velocity c(t), we must alter the 
Equivalency Principle and Field Equation so that locally there is a gener- 
alized Minkowski metric and a generalized Lorentz transform. We propose 
to alter the field equations by allowing c, G, and A to be functions of the 
cosmic time t and by introducing a generalized time t for which dt = c(t)dt, 
which yields the required generalization. This alteration shows that the 
equation for determining a(t) with a variable c(t) is the same as a(t) with 
c = 1. We have demonstrated that in the limit of a zero density universe 
the physical origin coordinates meet our physicality criteria for all distances 
with a constant photon velocity and become the well known coordinates first 
transformed from RW by Robinson [14]. 

The kinematics on distant galactic material should be calculated using 
a co- moving generalized Minkowski metric, whose photon velocity is c(t). 
Measurements of the fine structure in the absorption spectra from distant 
quasars caused by intervening clouds of matter show a change from that 
measured on earth. These measurements experimentally confirm that phys- 
ical "constants" do change with cosmic time. The measurements show a 
very small fractional change, far too small to be caused solely by the c(t) 
calculated in this paper. This implies that if we have a variable photon veloc- 
ity, some other physical "constants" should vary as well. We have followed 
Magueijo's approach[7] for a possible dependence of "constants" on c(t) in 
Sect 10, which indicate that the formulas for red shift and distance should 
be modified. The proposed modification explains the apparent universe ac- 
celeration found with supernovae la on high red-shift galaxies by assuming 
c(t) = vo(t) in a flat universe with no vacuum energy, and with the experi- 
mentally determined requirement that the SNIa luminosity be proportional 
to l/c(t). However, it does not eliminate the horizon problem. 

The agreement with the measurements of the Supernovae la suggest that 
these variable "constants" also be applied to other astrophysical observations 
where vacuum energy has been previously invoked, e.g., CMB, gravitational 
lensing, and dynamical estimates of cluster masses. It would also be a test of 
the theory if the secondary effect of the fine structure spectrum in the heavy 
metal atoms at high z could be calculated and compared to observation. 
Another challenge is to calculate the age of old globular clusters with c(t) 
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to see whether it is shortened enough to be still significantly less than the 
16% shorter age of a critical universe, 1/2Hq instead of 2/3Hq. 

This variation of c(t) should in principle be measurable by earth exper- 
iments. For a critical universe density of matter this fractional decrease 
should be about 1CT 14 in 10 4 sees, which should be measurable if the pro- 
jected fractional frequency accuracy of 10~ 17 in 10 4 sees is achieved, and 
if the frequency is proportional to c(t) as proposed in section 10. If the 
energy density due to the cosmological constant is dominant, the photon ve- 
locity and frequency will increase with time instead of decreasing, providing 
another test of the existence or non-existence of dark energy. 
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